Goulphar: rapid access and expertise for standard two-color microarray normalization methods

BACKGROUND: Raw data normalization is a critical step in microarray data analysis because it directly affects data interpretation. Most of the normalization methods currently used are included in the R/BioConductor packages but it is often difficult to identify the most appropriate method. Furthermore, the use of R commands for functions and graphics can introduce mistakes that are difficult to trace. We present here a script written in R that provides a flexible means of access to and monitoring of data normalization for two-color microarrays. This script combines the power of BioConductor and R analysis functions and reduces the amount of R programming required. RESULTS: Goulphar was developed in and runs using the R language and environment. It combines and extends functions found in BioConductor packages (limma and marray) to correct for dye biases and spatial artifacts. Goulphar provides a wide range of optional and customizable filters for excluding incorrect signals during the pre-processing step. It displays informative output plots, enabling the user to monitor the normalization process, and helps adapt the normalization method appropriately to the data. All these analyses and graphical outputs are presented in a single PDF report. CONCLUSION: Goulphar provides simple, rapid access to the power of the R/BioConductor statistical analysis packages, with precise control and visualization of the results obtained. Complete documentation, examples and online forms for setting script parameters are available from

Eoulsan : analyse du séquençage à haut débit dans le cloud et sur la grille

Eoulsan : analyse du séquençage à haut débit dans le cloud et sur la grill

Identification of Phox2b-regulated genes by expression profiling of cranial motoneuron precursors

<p>Abstract</p> <p>Background</p> <p>Branchiomotor neurons comprise an important class of cranial motor neurons that innervate the branchial-arch-derived muscles of the face, jaw and neck. They arise in the ventralmost progenitor domain of the rhombencephalon characterized by expression of the homeodomain transcription factors Nkx2.2 and Phox2b. Phox2b in particular plays a key role in the specification of branchiomotor neurons. In its absence, generic neuronal differentiation is defective in the progenitor domain and no branchiomotor neurons are produced. Conversely, ectopic expression of Phox2b in spinal regions of the neural tube promotes cell cycle exit and neuronal differentiation and, at the same time, induces genes and an axonal phenotype characteristic for branchiomotor neurons. How Phox2b exerts its pleiotropic functions, both as a proneural gene and a neuronal subtype determinant, has remained unknown.</p> <p>Results</p> <p>To gain further insights into the genetic program downstream of Phox2b, we searched for novel Phox2b-regulated genes by cDNA microarray analysis of facial branchiomotor neuron precursors from heterozygous and homozygous <it>Phox2b </it>mutant embryos. We selected for functional studies the genes encoding the axonal growth promoter Gap43, the Wnt antagonist Sfrp1 and the transcriptional regulator Sox13, which were not previously suspected to play roles downstream of <it>Phox2b </it>and whose expression was affected by <it>Phox2b </it>misexpression in the spinal cord. While <it>Gap43 </it>did not produce an obvious phenotype when overexpressed in the neural tube, <it>Sfrp1 </it>induced the interneuron marker Lhx1,5 and <it>Sox13 </it>inhibited neuronal differentiation. We then tested whether <it>Sfrp1 </it>and <it>Sox13</it>, which are down-regulated by Phox2b in the facial neuron precursors, would antagonize some aspects of <it>Phox2b </it>activity. Co-expression of <it>Sfrp1 </it>prevented <it>Phox2b </it>from repressing Lhx1,5 and alleviated the commissural axonal phenotype. When expressed together with <it>Sox13</it>, <it>Phox2b </it>was still able to promote cell cycle exit and neuronal differentiation, but the cells failed to relocate to the mantle layer and to extinguish the neural stem cell marker Sox2.</p> <p>Conclusion</p> <p>Our results suggest novel roles for <it>Sfrp1 </it>and <it>Sox13 </it>in neuronal subtype specification and generic neuronal differentiation, respectively, and indicate that down-regulation of <it>Sfrp1 </it>and <it>Sox13 </it>are essential aspects of the genetic program controlled by Phox2b in cranial motoneurons.</p

Kinetic transcriptome analysis reveals an essentially intact induction system in a cellulase hyper-producer Trichoderma reesei strain

International audienceBackground: The filamentous fungus Trichoderma reesei is the main industrial cellulolytic enzyme producer. Several strains have been developed in the past using random mutagenesis, and despite impressive performance enhancements, the pressure for low-cost cellulases has stimulated continuous research in the field. In this context, comparative study of the lower and higher producer strains obtained through random mutagenesis using systems biology tools (genome and transcriptome sequencing) can shed light on the mechanisms of cellulase production and help identify genes linked to performance. Previously, our group published comparative genome sequencing of the lower and higher producer strains NG 14 and RUT C30. In this follow-up work, we examine how these mutations affect phenotype as regards the transcriptome and cultivation behaviour. Results: We performed kinetic transcriptome analysis of the NG 14 and RUT C30 strains of early enzyme production induced by lactose using bioreactor cultivations close to an industrial cultivation regime. RUT C30 exhibited both earlier onset of protein production (3 h) and higher steady-state productivity. A rather small number of genes compared to previous studies were regulated (568), most of them being specific to the NG 14 strain (319). Clustering analysis highlighted similar behaviour for some functional categories and allowed us to distinguish between induction-related genes and productivity-related genes. Cross-comparison of our transcriptome data with previously identified mutations revealed that most genes from our dataset have not been mutated. Interestingly, the few mutated genes belong to the same clusters, suggesting that these clusters contain genes playing a role in strain performance. Conclusions: This is the first kinetic analysis of a transcriptomic study carried out under conditions approaching industrial ones with two related strains of T. reesei showing distinctive cultivation behaviour. Our study sheds some light on some of the events occurring in these strains following induction by lactose. The fact that few regulated genes have been affected by mutagenesis suggests that the induction mechanism is essentially intact compared to that for the wild-type isolate QM6a and might be engineered for further improvement of T. reesei. Genes from two specific clusters might be potential targets for such genetic engineering

Insights into the role of DNA methylation in diatoms by genome-wide profiling in Phaeodactylum tricornutum

DNA cytosine methylation is a widely conserved epigenetic mark in eukaryotes that appears to have critical roles in the regulation of genome structure and transcription. Genome-wide methylation maps have so far only been established from the supergroups Archaeplastida and Unikont. Here we report the first whole-genome methylome from a stramenopile, the marine model diatom Phaeodactylum tricornutum. Around 6% of the genome is intermittently methylated in a mosaic pattern. We find extensive methylation in transposable elements. We also detect methylation in over 320 genes. Extensive gene methylation correlates strongly with transcriptional silencing and differential expression under specific conditions. By contrast, we find that genes with partial methylation tend to be constitutively expressed. These patterns contrast with those found previously in other eukaryotes. By going beyond plants, animals and fungi, this stramenopile methylome adds significantly to our understanding of the evolution of DNA methylation in eukaryotes.Fil: Veluchamy, Alaguraj. Institut de Biologie de l'École Normale Supérieure; FranciaFil: Lin, Xin. Institut de Biologie de l'École Normale Supérieure; Francia. Xiamen University; ChinaFil: Maumus, Florian.Fil: Rivarola, Maximo Lisandro.Fil: Bhavsar, Jaysheel.Fil: Creasy, Todd.Fil: O'Brien, Kimberly.Fil: Sengamalay, Naomi A..Fil: Tallon, Luke J..Fil: Smith, Andrew D..Fil: Rayko, Edda.Fil: Ahmed, Ikhlak.Fil: Crom, Stéphane Le.Fil: Farrant, Gregory K..Fil: Sgro, Jean-Yves.Fil: Olson, Sue A..Fil: Bondurant, Sandra Splinter.Fil: Allen, Andrew.Fil: Rabinowicz, Pablo D..Fil: Sussman, Michael R..Fil: Bowler, Chris.Fil: Tirichine, Leïla

Erratum to: Genome sequencing of the Trichoderma reesei QM9136 mutant identifies a truncation of the transcriptional regulator XYR1 as the cause for its cellulase-negative phenotype

International audienceBackground : Trichoderma reesei is the main industrial source of cellulases and hemicellulases required for the hydrolysis of biomass to simple sugars, which can then be used in the production of biofuels and biorefineries. The highly productive strains in use today were generated by classical mutagenesis. As byproducts of this procedure, mutants were generated that turned out to be unable to produce cellulases. In order to identify the mutations responsible for this inability, we sequenced the genome of one of these strains, QM9136, and compared it to that of its progenitor T. reesei QM6a.Results : In QM9136, we detected a surprisingly low number of mutagenic events in the promoter and coding regions of genes, i.e. only eight indels and six single nucleotide variants. One of these indels led to a frame-shift in the Zn2Cys6 transcription factor XYR1, the general regulator of cellulase and xylanase expression, and resulted in its C-terminal truncation by 140 amino acids. Retransformation of strain QM9136 with the wild-type xyr1 allele fully recovered the ability to produce cellulases, and is thus the reason for the cellulase-negative phenotype. Introduction of an engineered xyr1 allele containing the truncating point mutation into the moderate producer T. reesei QM9414 rendered this strain also cellulase-negative. The correspondingly truncated XYR1 protein was still able to enter the nucleus, but failed to be expressed over the basal constitutive level.Conclusion : The missing 140 C-terminal amino acids of XYR1 are therefore responsible for its previously observed auto-regulation which is essential for cellulases to be expressed. Our data present a working example of the use of genome sequencing leading to a functional explanation of the QM9136 cellulase-negative phenotype

An evaluation of custom microarray applications: the oligonucleotide design challenge

The increase in feature resolution and the availability of multipack formats from microarray providers has opened the way to various custom genomic applications. However, oligonucleotide design and selection remains a bottleneck of the microarray workflow. Several tools are available to perform this work, and choosing the best one is not an easy task, nor are the choices obvious. Here we review the oligonucleotide design field to help users make their choice. We have first performed a comparative evaluation of the available solutions based on a set of criteria including: ease of installation, user-friendly access, the number of parameters and settings available. In a second step, we chose to submit two real cases to a selection of programs. Finally, we used a set of tests for the in silico benchmark of the oligo sets obtained from each type of software. We show that the design software must be selected according to the goal of the scientist, depending on factors such as the organism used, the number of probes required and their localization on the target sequence. The present work provides keys to the choice of the most relevant software, according to the various parameters we tested

The CRE1 carbon catabolite repressor of the fungus Trichoderma reesei: a master regulator of carbon assimilation

<p>Abstract</p> <p>Background</p> <p>The identification and characterization of the transcriptional regulatory networks governing the physiology and adaptation of microbial cells is a key step in understanding their behaviour. One such wide-domain regulatory circuit, essential to all cells, is carbon catabolite repression (CCR): it allows the cell to prefer some carbon sources, whose assimilation is of high nutritional value, over less profitable ones. In lower multicellular fungi, the C2H2 zinc finger CreA/CRE1 protein has been shown to act as the transcriptional repressor in this process. However, the complete list of its gene targets is not known.</p> <p>Results</p> <p>Here, we deciphered the CRE1 regulatory range in the model cellulose and hemicellulose-degrading fungus <it>Trichoderma reesei </it>(anamorph of <it>Hypocrea jecorina</it>) by profiling transcription in a wild-type and a delta-<it>cre1 </it>mutant strain on glucose at constant growth rates known to repress and de-repress CCR-affected genes. Analysis of genome-wide microarrays reveals 2.8% of transcripts whose expression was regulated in at least one of the four experimental conditions: 47.3% of which were repressed by CRE1, whereas 29.0% were actually induced by CRE1, and 17.2% only affected by the growth rate but CRE1 independent. Among CRE1 repressed transcripts, genes encoding unknown proteins and transport proteins were overrepresented. In addition, we found CRE1-repression of nitrogenous substances uptake, components of chromatin remodeling and the transcriptional mediator complex, as well as developmental processes.</p> <p>Conclusions</p> <p>Our study provides the first global insight into the molecular physiological response of a multicellular fungus to carbon catabolite regulation and identifies several not yet known targets in a growth-controlled environment.</p

Genome-Wide Transcriptome Analyses of Silicon Metabolism in Phaeodactylum tricornutum Reveal the Multilevel Regulation of Silicic Acid Transporters

BACKGROUND:Diatoms are largely responsible for production of biogenic silica in the global ocean. However, in surface seawater, Si(OH)(4) can be a major limiting factor for diatom productivity. Analyzing at the global scale the genes networks involved in Si transport and metabolism is critical in order to elucidate Si biomineralization, and to understand diatoms contribution to biogeochemical cycles. METHODOLOGY/PRINCIPAL FINDINGS:Using whole genome expression analyses we evaluated the transcriptional response to Si availability for the model species Phaeodactylum tricornutum. Among the differentially regulated genes we found genes involved in glutamine-nitrogen pathways, encoding putative extracellular matrix components, or involved in iron regulation. Some of these compounds may be good candidates for intracellular intermediates involved in silicic acid storage and/or intracellular transport, which are very important processes that remain mysterious in diatoms. Expression analyses and localization studies gave the first picture of the spatial distribution of a silicic acid transporter in a diatom model species, and support the existence of transcriptional and post-transcriptional regulations. CONCLUSIONS/SIGNIFICANCE:Our global analyses revealed that about one fourth of the differentially expressed genes are organized in clusters, underlying a possible evolution of P. tricornutum genome, and perhaps other pennate diatoms, toward a better optimization of its response to variable environmental stimuli. High fitness and adaptation of diatoms to various Si levels in marine environments might arise in part by global regulations from gene (expression level) to genomic (organization in clusters, dosage compensation by gene duplication), and by post-transcriptional regulation and spatial distribution of SIT proteins